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ABSTRACT 

Dielectronic recombination (DR) of singly charged ions is a reaction pathway 
that is commonly neglected in chemical models of molecular clouds. In this 
study we include state-of-the-art DR data for He+, C"*", N+, O"*", Na"*", and Mg"*" 
in chemical models used to simulate dense molecular clouds, protostars, and 
diffuse molecular clouds. We also update the radiative recombination (RR) rate 
coefficients for H+, He^, C^, N"^, 0+, Na^, and Mg+ to the current state-of- 
the-art values. The new RR data has little effect on the models. However, the 
inclusion of DR results in significant differences in gas-grain models of dense, cold 
molecular clouds for the evolution of a number of surface and gas-phase species. 
We find differences of a factor of 2 in the abundance for 74 of the 655 species at 
times of 10^-10^ years in this model when we include DR. Of these 74 species, 
16 have at least a factor of 10 difference in abundance. We find the largest 
differences for species formed on the surface of dust grains. These differences are 
due primarily to the addition of C"*" DR, which increases the neutral C abundance, 
thereby enhancing the accretion of C onto dust. These results may be important 
for the warm-up phase of molecular clouds when surface species are desorbed into 
the gas phase. We also note that no reliable state-of-the-art RR or DR data exist 
for Si+, P"*", S'^, Cl"^, and Fe+. Modern calculations for these ions are needed to 
better constrain molecular cloud models. 
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Introduction 



Investigating the physics and chemistry of molecular clouds (also known as quiescent 
cores) is crucial if one is to understand the processes that ultimately lead to star forma- 
tion. Such studies are also important in the field of astrobiology since these clouds are the 
birthplace of the first organic molecules. 

Chemical models used to describe the evolution of molecular clouds typically include 
hundreds of species. At molecular cloud temperatures (< 100 K) these atoms and molecules 
are either neutral or singly ionized. The abundances of these species are, in turn, governed 
by thousands of reactions that are highly non-linear. These derived abundances are sen- 
sitive to the accuracy of the rate coefficients involved. The implicatio n of uncertai i ities i n 
the rate coefficients used iii chemical models has been in vestigated by iRoueff et al.l (Il996l ) , 
Vasvunin et"aD J2004j . boosh . IWakelam et al.l J2005I . boooh . and others. However, these stud- 
ies cannot account for reaction processes that are not included in the models in the first 
place. 

To the best of our knowledge, the only recombination process of free electrons with 
atomic ions included in any rnolecul ar cloud simulation is radiative recombination (RR). 
Recent calculations by Badnell ( 2006 1^ have improved the accuracy of the rate coefficients of 
these reactions for a number of ions. The alternative pathway of dielectronic recombination 
(DR) has previously not been included in molecular cloud models. This is largely due 
to the bulk of published DR calculations and experiments being v alid only for p lasma s 
of higher temperature. Recently, however, theoretical calculations by Badnell et al. (j2003 l^ 
have probed DR at the low temperature regimes of molecular clouds. At these temperatures, 
the DR rate coefficient can be a factor of 5 or more greater than the RR rate coefficient for 
certain systems. 

In the present paper, we both update the RR data and include DR for the relevant singly 
charged atomic ions. We do this for several molecu lar cloud models under a variety of initial 
conditions. We use Naho on (IWakelam et al.ll2004r) to simulate derise clo uds and protostars, 
the OSU gas-grain code (iHasegawa et al.l I1992I : I Garrod fc HerbstI 120061) to sim ulate dense 
clouds, and the Meudon Photodissociation Region (PDR) code ( lie Petit et al.l |2002| . |2006| ) 
to simulate both diffuse and dense PDRs. We compare the species abundances with those 
calculated with and without DR included in order to show the effects on the models. The 
remainder of the paper is organized as follows: In Sec. [2] we review the recent developments 



^http: / / amdpp. phys.strath.ac.uk/tamoc/RR/ 
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in our understanding RR and DR. Section [3] outlines the chemical models we use to simulate 
dense clouds, diffuse clouds, and protostars. In Sec. H] we present the results of including DR 
in these chemical models. Concluding remarks are given in Sec. O 



Recombination data 



2.1. Radiative Recombination 



Radiative recombination (RR) is a one-step recombination process that occurs when a 
free electron is captured by an ion. Energy and momentum are conserved in the process 
by the simultaneous emission of a photon. At the low temperatures typical of molecular 
clouds, RR is the dominant ion-electron recombination process for most ions. Until now, 
the molecular cloud codes that we use in the present work (see Sec. |3D have used RR rate 
coefficients from the UMIST databas^ ( Woodall et al. 2007 ). For the most part it is unclear 
where the RR data in the UMIST database stem from as there is often no reference given; this 
is the case for the rate coefficients of H"*", He"'", Na"*", Mg"*", Si^, P"*", S"*", Cl "^, and Fe"*". There 
is no RR data for F+. For C"*" and the rate coefficients are those of Nahar fc Pradhan 



( 119971 ) and for O"*" the rate coefficient is from lNahar fc PradhanI (119991 ). The data of Nahar & 
Pradhan are actually unified RR-I-DR calculations using LS-coupling. As a result, they have 
no DR component at the low temperatures of molecular clouds since they do not account 
for fine structure transitions of the ground term (see Sec. 12.21) . 

In recent years there have been atte mpts to better understand the RR process using 
state-of-the-art computational techniques. iBadnelll (120061 ) has calculated RR rate coefficients 
for all elements from H through Zn for bare through Mg-like isoelectronic sequences. We have 
implemented these calculations for the ions included in the chemical networks of the models 
consid ered here, namely II+, IIe+, C+, N+, O"^, Na^, and Mg"*". For the ions that iBadnell 
( 120061 ) has not ca lculated — namel y Si+, S+, CI"*", and Fe"*" — we have used the rate coefficients 



recommended by iMazzotta et al.l (119981 ). These come from lAldrovandi fc Pequignotl (119731 ) 
for Si"*" and S"*"; an extrapolation of the S-like sequence is used to estimate the CI"*" rate 
coefficient; and the rate coefficient from Schull & van Steenberg (1982; refitted with the 
Verner & Ferland 1996 formula) is used for Fe+. 

At molecular cloud temperatures, th e largest differ ence found between the RR rate 
coefficients of the UMIST dataset and the iBadnelll (120061 ) RR data is for Mg+, as is shown 
in Fig. m At a temperature of 10 K, typical of a cold molecular cloud, the UMIST rate 
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coefficient is 60% larger than that of Badnell. 



2.2. Dielectronic recombination 

Dielectronic recombination (DR) is a two-step recombination process that begins when 
a free electron coUisionally excites a core electron of an ion and is simultaneously captured. 
The core electron excitation can be labeled nlj '^'^'j'l where n is the principal quantum 
number, / the orbital angular momentum, and j the total angular momentum. The energy 
of this intermediate system lies in the continuum and the complex may autoionize. The 
DR process is complete when the system emits a photon, reducing the total energy of the 
recombined system to below its ionization threshold. Conservation of energy requires that 
for DR to go forward 

Ek = AE- Eb. (1) 

Here E^ is the kinetic energy of the incident electron, AE the excitation energy of the 
initially bound electron in the presence of the captured electron, and Ej, the binding energy 
released when the incident electron is captured onto the excited ion. Because AE and Ei, 
are quantized, DR is a resonant process occurring for a given channel at energies Ek < AE. 

For a molecular cloud of temperature Tdoud, the important DR channels are those where 
AE ~ ksTcioud ^ 0.01 eV. Thus, from atomic energy structure considerations alone, it is clear 
that fine structure excitations of the ground term are important in molecular clouds. Until 
recently, the preponderance of DR data have been calculated using LS-coupling. However, 
this coupling s cheme does no t inclu de these low-energy resonances. Recent state-of-the-art 



calculations by iBadnell et al.l (120031 ) . using intermediate coupling, have accounted for these 
fine structure channels, calculating DR rate coefficients for H- through Mg-like ions of all 
elements from He through Zn. 

For the singly-ionized ions of interest in molecular clouds — H"*", He+, C+, N+, 0+, F+, 
Na+, Mg+, Si"*", P"*", S+, C1+, and Fe+ — one would expect a significant low-temperature DR 
contribution for C"*", N"*", F"*", Si^, P^, Cl^ and Fe"*" due to their fine structure splitting of 
the ground term. The other ions either do not undergo DR (i.e., H"*") or have ii o such fine 



structure splitting. Of these five ions with fine structure splitting, iBadnell et al.l (120031 ) has 
calculated DR only for C"*", N+, and F+. Figs. [2] and [3] show a comparison of the DR and RR 
rate coefficients for C"^ and N"*", respectively. At temperatures < 100 K, typical of molecular 
clouds, the DR component dominates the electron-ion recombination rate coefficient high- 
lighting the importance of including this reco mbination pro c ess. F or the present work we 



have implemented the DR rate coefficients of IBadnell et al.l (120031 ) for He+, C"*", N"*", O 



F+, Na^, and Mg"*". No such non-LS coupling calculations exist for Si"*", P'^, Cl^ and Fe^ 
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Models 



3.1. Nahoon 



The Nahoon code (jWakelam et al.l |2004| ) is a pseudo-time-dependent chemical model 
that computes the chemical evolution of gas-phase species for a fixed gas temperature and 
density. It is suited for modeling dense molecular clouds and protostars. The model includes 
452 species with initial conditions being atomic (neutral and singly charged) He, C, N, O, 
F, Na, Mg, Si, P, S, CI, and Fe plus molecular hydrogen. Dust grains are included in 
the model but are not fully treated; ions are allowed to neutralize by charge exchange with 
negative grains but accretion of species onto grains is not included, and thus no grain surface 
chemistry is possible. The code computes the evolution of the abundance of each species 
as governed by the 4423 included reactions between species. T he reaction rate coefficients 
come from the osu. 20050 chemical network ( ISmith et al.ll2004j ). The RR rate coefficients 
in this network are those from the UMIST data base (see discu ssion in Sec. 12. ip . We have 
replaced these RR rate coefficients with those of iBadnelll (120061 ) for H+, He+, N^, O"*", 
Na+, and Mg+. We have also added the PR react ions for He+, C+, N+, 0+, Na+, and Mg+ 
using the rat e coefficients of iBadnell et al.l (|2003r). For the other ions that have not been 
calculated bv lSadnell eTaP J2003h orlSadnelll J2006 h — Si+ , S+, C1+, and Fe+— we have used 



the RR and DR rate coefficients of Mazzotta et al 



(119981 ). While the elements F and P are 



present in Nahoon, recombination of F+ and P+ are not included. We have chosen not to 
introduce these processes and instead concentrate on the effect of the updated RR and DR 
rate coefficients for those reactions that are already present. 

With these new reactions included, we have used Nahoon to simulate the evolution of 
dense cold clouds and protostars. A dense cold cloud is simulated by assuming a temperature 
of 10 K, an initial H2 density of 10^ cm~^, a visual extinction of 10 so that the photochemistry 
driven by external UV photons does not occur, and a fixed cosmic-ray ionization rate of 1.3 x 
10"^'' s~^. Th ese dense cloud condit ions are given in Tabled! We use the low-metal elemental 
abundances of iGraedel et al.l (119821 ) as our initial conditions. These are also given in Tabled] 
For a protostar, we take our initial gas-phase chemical composition as that computed by 
Nahoon for the above dense cloud conditions at 10^ years. Additionally, we increased the 
temperature to 100 K, the initial H2 density to 10^ cm~^, and set the abundances of H2O, 
H2CO, CH4O, and CH4 as given in Table dl This is required because Nahoon does not 
include the surface chemistry that is needed to form these species, so we artificially increase 
their abundance in the gas phase to simulate protostar conditions. Results for both dense 



http://www.physics.ohio-state.edu/~eric/researchJiles/osu. 2005 
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clouds and protostars are presented in Sec. |H 



3.2. OSU gas-grain code 



The OSU gas-grain code (jHasegawa et al.l Il992l : iGarrod fc HerbstI l2006l ) is similar to 
Nahoon in that it follows the evolution of the species abundances as a function of time for 
a fixed cloud temperature and density. Like Nahoon, it is used to simulate the conditions of 
dense molecular clouds. However, unlike Nahoon, the OSU gas-grain code includes surface 
chemistry in addition to gas-phase chemistry. This allows for time-dependent accretion onto 
grain surfaces as well as thermal and cosmic ray-induced evaporation from the surfaces. The 
code assumes a sticking coefficient of unity for neutral species that s trike grains. It is unclear 
whether charged species can also accrete onto grains (IWatsonlll976l ). It is assumed here that 
they do not, with the energy released on neutralization resulting in the species desorbing from 
the grain surface. T his model also includes the new non-thermal evaporation mechanism of 
Garrod et al.l (120071 ). This process assumes that the energy released by exothermic surface 
reactions partially evaporate the products. 

Species can react on the grain surfaces to form molecules more efficiently than is possible 
solely in the gas phase. This surface chemistry is an important contributor to the species 
present in the warm-up phase of hot molecular cores (IGarrod &: HerbstI l2006l ). Species that 
form on the grains at earlier times are injected into the gas phase when the dust is heated 
by a nearby protostar. 

The OSU code uses the osu.2005 chemical network, although there is no F present in 
the OSU gas-grain code, unlike in Nahoon. With the addition of the surface species to the 
network used with Nahoon, there are 655 species and 6309 reactions. We added the new RR 
and DR rate coefficients in the same way as was done in Nahoon (see Sec. 13. ip . With these 
changes to the chemical network, we have run the OSU code with the same initial conditions 
as used in Sec. 13. II when simulating a dense cold cloud (see Table [T]). We present our results 
in Sec. ID 



3.3. PDR code 



In order to simulate a molecular clou d with a significant impin gent radiation field, we 
use the PDR code of the Meudon group (lie Petit et al.ll2002l . I2006I 1H This code is used to 



^http://aristote. obspm.fr/MIS/pdr/pdrl. html 
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simulate clouds with no rapidly evolving processes occurring. The cloud is treated in steady 
state as an infinite slab of gas and dust irradiated by an ultraviolet radiation field impinging 
on both sides of the cloud. While dust is included in the model, there is no surface chemistry 
other than allowing for the formation of H2 which cannot be produced rapidly enough in the 
gas phase. 

Using a standard chemistry file with 120 species we have added the new RR and DR 
data to the PDR code. The elemental abundances are given in Table [H We have run the 
model with physical conditions that simulate a diffuse PDR and, separately, a dense PDR. 
For diffuse conditions, we set the H2 density to 25 cm~^, the temperature to 20 K, and the 
cosmic ray ionization rate of 5 x 10~^^ s~^. The radiation field is 1 Draine at the edges 



of the cloud. For a dense PDR, we simulate the Horsehead nebula (jPety et al.ll2005l ) with 
an H2 density of 5 x 10^ cm~^, a temperature of 90 K, and a cosmic ray ionization rate 
of 5 X 10"^^ s"^ For this case the radiation field is 100 Draine at the edge of the cloud. 
The species abundances are calculated as a function of depth into the cloud. Results are 
presented in Sec. HI 



4. Results 

For each of the models described in Sec. Owe have run the simulation with the original 
chemical network without any additions, with just the new RR, and again with the new RR 
and DR rate coefficients included. When running the models with only the new RR data 
included we find no differences greater than a factor of 2 for any of the models. Only for 
the OSU gas-grain code do we find abundance differences of greater than 50%. These are 
found for gas-phase Fe, surface FeH, and surface MgH2, resulting from the change in the 
Mg"*" and Fe"*" RR rate coefficients. The more significant results reported in the remainder 
of this section can thus be attributed primarily to the effect of DR on the models. 

Simulations of a dense, cold molecular cloud of low metallicity were run with the OSU 
gas-grain code, which computes both gas-phase and surface chemistry. We find that the 
introduction of DR results in significant differences in certain species evolution. To identify 
those differences that are most important in the cloud chemistry we consider differences only 
at times when species have an abundance (either before or after the inclusion of DR) of at 
least 10^^^ with respect to the total H nuclei density. This represents those species that are 
most likely to be detectable. Secondly, we only consider abundances at times of 10^-10^ years 
in the model evolution as it is during this period in the evolution of the model that observed 
molecular cloud conditions are best represented. With these criteria, we find 100 species to 
have a difference in abundance of at least 50% on including DR. We list in Table [2] the 74 
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species for which we find at least a factor of 2 difference in abundance. Of these 74 species 
48 are attached to the surface of dust grains and 26 are found in the gas phase. The effect 
of the new DR rate coefficients on species abundances is over 2 orders of magnitude for 16 
species. The greatest effect is found for surface O3, with a 3 orders of magnitude difference 
in abundance found. In Figs. HJ [5l and [6] we show the evolution of the abundance of surface 
C, surface CH3OH, and surface O3, respectively. 

We have run the gas-grain model a number of times, including DR for only a single 
selected ion in each run, and repeated this for every ion. This allows us to determine which 
ion has the largest effect on the differences we have found. These results indicate that the 
abundance differences listed in Table [2] are primarily due to the introduction of the C"*" DR. 
When the C"*" DR rate coefficient is not included in the model we find no differences in 
abundance greater than a factor of 2. The inclusion of DR increases the C"*" recombination 
rate, neutralizing C earlier in the model evolution. The model allows neutral species to 
accrete onto grains but not charged species. Thus, with an increased abundance of neutral 
C, the surface chemistry is enriched and results in the abundance differences that we observe. 

Despite the large increase in the electron-ion recombination rate coefficient for N"*" 
(Fig. [3]), the comparatively small effect of the inclusion of N"*" DR can be explained by 
considering the ionization balance of N in the simulation. Initially (see Table [1]) all N in 
the model is in the neutral charge state, so recombination of N"*" with a free electron is not 
possible. As the model evolves, some N"*" is formed but neutral N remains many orders of 
magnitude greater in abundance. 

In contrast, atomic C, due to its lower first ionization potential, is initially all singly 
ionized by UV radiation shortward of the 13.6 eV ionization potential of H but above the 
11.3 eV ionization potential of C (again see Tabled]). Thus, recombination of C"*" is important 
from the very start of the model evolution and the inclusion of DR greatly increases the 
electron-ion recombination rate. The increase in the C recombination rate also results in a 
change in the ionization fraction of the gas as can be seen in the abundance evolution of the 
sum of all negatively charged species. To show this we plot the free electron abundance in 
Fig. [71 The abundance of all other negatively charged species are negligible in comparison. 

To investigate the importance of the surface chemistry on our results, we have run the 
gas-grain model with no accretion of species onto dust grains allowed. We have compared 
results with and without DR included. With the surface chemistry excluded, we find differ- 
ences in abundance of 50% for 32 species but factor of 2 increases for only 5 species — CH3CN, 
CH3CO"'", HC2NC, C3HJI", and HC3N. These 5 gas-phase molecules had among the smallest 
changes of the 74 species we identified as having a factor of 2 difference when the surface 
chemistry was included. Of the 74 species listed in Table [2l 48 are surface species which are 
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unable to form when we do not include the surface chemistry. However, there are also 26 
gas-phase species in Table [21 Only the 5 gas-phase species identified above have differences 
of a factor of 2 when the surface chemistry is excluded. Of the 21 other gas-phase species in 
Tabled 11 have no formation routes without the surface chemistry included. The other 10, 
whose abundances are greatly enhanced by formation routes on grain surfaces, have changes 
of less than 50% in a purely gas-phase model. Thus, the majority, and indeed the largest, 
abundance differences seen when including DR in the OSU gas-grain model are found for 
species influenced by surface reactions. 

Not surprisingly then, we do not find significant abundance differences when running 
the dense cloud simulation using Nahoon, which does not include grain surface chemistry. 
All species in Nahoon are in the gas phase. The introduction of the DR data resulted in 20 
species having an abundance difference of 50% and only 2 species, OCS and C2H3N, having 
abundance increases of greater than a factor of 2. When running the protostar simulation 
using Nahoon we again find few significant differences when DR is included; 21 species have 
an abundance difference of 50% or more. The only species with greater than a factor of 2 
difference are C2S, HON, HNC, C2H3N, HCOOCH3, CH3OCH3, H2CN+, and CH3OCHI, 
all of which are increased in abundance with DR relative to without. 

Simulations of diffuse and dense PDRs were carried out using the Meudon PDR code, 
allowing us to probe the effect of DR in a molecular cloud with a radiation field present. On 
introducing the new DR rate coefficients there is a large increase in the total recombination 
rate of C"*" and we see a significant effect on the abundance of neutral carbon as a result. 
This is the case for both diffuse and dense PDRs. The column density of C is shown in 
Figs. [8] and [9] for diffuse and dense regions, respectively, as a function of visual extinction 
(i.e., depth into the cloud). However, we find no significant abundance differences for any 
other species in the diffuse case. For a dense PDR, there are only 5 other species that show 
an abundance difference greater than 50% — SO2, O2H+, HSO2 , HOCS+, and C4. These 
species all increase in abundance with the introduction of DR but the increases are all less 
than a factor of 2. 



5. Conclusions 

This work has investigated the importance of new RR and DR rate coefficients in molec- 
ular cloud models. We have updated the RR rate coefficients of H+, He"*", C"*", N+, O'^, Na^, 
and Mg"^ to the current state-of-the-art in 3 different molecular cloud chemical models. We 
have run these models under different conditions and compared the abundances of the species 
present with and without the new RR rate coefficients. We find that the new RR data have 
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no significant effects on tfie model results (i.e., less than a factor of 2 effect). We have also 
included DR for He"*", C"*", N"*", O"*", Na"*", and Mg"*" in the 3 chemical models, finding some 
sizable differences (greater than a factor of 2) in the evolution of certain species. 

We find that DR has the greatest effect for dense cloud models, particularly when surface 
chemistry is included. Of the models we ran here, only the OSU gas-grain code includes the 
extensive surface chemistry that allows species to be accreted onto grains and subsequently 
react with other surface species. Nahoon does not include this surface chemistry and the 
Meudon code only allows the formation of H2 on surfaces. DR has its greatest influence 
for species that form on grain surfaces. These species may be particularly important in the 
warm-up phase of molecular clouds when they become dcsorbed into the gas phase. We 
find DR of C"^ to be primarily responsible for the abundance changes in the model. This is 
because charged species are assumed not to accrete onto grains and adding DR increases the 
gas-phase neutral C abundance available for first accreting onto grains and then participating 
in the surface chemistry. 

We conclude by noting that there are some ions for which no state-of-the-art RR or DR 
data exist, specifically Si+, P+, S"^, C1+, and Fe"*". Given the fine structure present in the 
ground terms of Si"^, Cl"^, P"*", and Fe"*" one would expect a significant DR component to the 
total electron-ion recombination rate coefficient at low temperatures for these ions. These 
4 elements also have first ionization potentials below that of H (13.6 cV). Given that they 
would thus be ionized in molecular clouds, the effect of increasing their recombination rates 
may be important. Modern calculations of both RR and DR rate coefficients for these ions 
are needed for generating reliable molecular cloud chemical models. 
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Table 1. Initial species abundances relative to the total H nuclei density, and physical 
conditions for dense clouds, protostars, and dense and diffuse PDRs. 



Mq T* Q TYl p+ pT* 
± Chi. dlllC' ud 




Initial conditions 






Protostar^ 


Diffuse PDR 




H 


— 


— 


9.90 X 10-^ 


9.90 X 10-1 


He 


1.40 X 10^^ 


1.40 X 10^^ 


1.00 X 10^1 


1.00 X 10-1 


N 


2.14 X 10"^ 


2.14 X 10-5 


7.50 X 10-5 


7.50 X 10-5 





1.76 X 10"^ 


1.76 X 10-^ 


3.19 X 10-^ 


3.19 X 10-^ 


F 


6.68 X 10-9 


6.68 X 10-9 


— 





C+ 


7.30 X 10-5 


7.30 X 10-5 


1.32 X 10-^ 


1.32 X 10-^ 


Na+ 


2.00 X lO-'' 


2.00 X 10-9 


— 


— 


Si+ 


8.00 X lO-'' 


8.00 X 10-9 




— 


s+ 


8.00 X 10-s 


O.UU X iU 


l.oO X lU 


1.86 X 10-5 


Mg+ 


7.00 X 10-9 


7.00 X 10-9 




- 


P+ 


3.00 X 10-9 


3.00 X 10-9 






C1+ 


4.00 X 10-9 


4.00 X 10-9 






Fe+ 


3.00 X 10-9 


3.00 X 10-9 


1.50 X 10-s 


1.50 X 10-^ 


H2O 




5.00 X 10-5 






H2CO 




2.00 X 10-6 






CH4O 




2.00 X 10-6 






CH4 




5.00 X 10-^ 






Temperature (K) 


10 


100 


20 


90 


H2 density (cm~^) 


1 X 10^ 


1 X 10^ 


25 


5 X 10^ 


Cosmic ray ionization rate 


1.3 X 10-^^ 


1.3 X 10-1^ 


5.0 X 10-^^ 


5.0 X 10-1^ 


Visual extinction 


10 


10 


1 


10 



^Thc elemental abundances arc those at the beginning of the simulation, while the 
molecular abundances (H2O, H2CO, CH4O, and CH4) are set to the given values at a 
model time of 105 years. 



Note. ~- 



There is no F included in the OSU gas-grain code. 



-13- 



Table 2. Maximum abundance differences found wlien including DR in tlie OSU gas-grain 

code relative to excluding DR. 







Maximum 


Time 


Relative 


Gas- Phase 


Surface 


Species 


Phase 


Difference 


(years) 


Abundance 


Detectioni 


Detection- 


O3 


Surface 


1440 


1.78 


X 


10^ 


1.06 X 10"^ 


No 


No 


N2H2 


Surface 


59.3 


3.16 


X 


10^ 


1.47 X 10"^^ 


No 


No 


H2O2 


Surface 


45.0 


5.62 


X 


10^ 


7.23 X 10 ^ 


No 


No 


HNCO 


Gas 


38.7 


1.78 


X 


10^ 


6.13 X 10"^^ 


Yes 


No 


HNCO 


Surface 


37.0 


1.78 


X 


10^ 


2.93 X 10 ^ 


Yes 


No 


NO 


Surface 


34.9 


3.16 


X 


10^ 


5.44 X 10~^^ 


Yes 


No 


NH2CHO 


Surface 


30.0 


5.62 


X 


10^ 


3.46 X 10 


Yes 


No 


C2H 


Surface 


27.2 


1.00 


X 


10^ 


1.06 X 10 


Yes 


No 


C2H4 


Surface 


23.1 


1.00 


X 


10^ 


5.64 X 10"^° 


Yes 


No 


HC2O 


Gas 


18.3 


1.78 


X 


10^ 


2.41 X 10"^^ 


No 


No 


/^x X 

CH2CO 


Surface 


15.0 


1.00 


X 


10^ 


3.10 X 10~^° 


No 


No 


CH2OH 


Gas 


14.1 


1.00 


X 


10^ 


5.47 X 10 


No 


No 


CH3OH 


Gas 


13.8 


1.00 


X 


10^ 


3.26 X 10 


Yes 


Yes 


C3H3N 


Surface 


12.8 


3.16 


X 


10^ 


1.98 X 10 


No 


No 


HC3N 


Surface 


12.5 


3.16 


X 


10^ 


3.77 X 10~^^ 


Yes 


No 


C2H6 


Gas 


11.3 


1.00 


X 


10^ 


5.07 X 10~^^ 


No 


No 


CH2NH2 


Gas 


9.76 


5.62 


X 


10^ 


1.96 X 10 


No 


No 


CH3NH 


Gas 


9.76 


5.62 


X 


10^ 


3.65 X 10 


No 


No 


CH5N 


Surface 


9.55 


5.62 


X 


10^ 


9.48 X 10"^ 


No 


No 


C2H2 


Surface 


8.27 


1.00 


X 


10^ 


1.19 X 10"^ 


Yes 


No 


CH3OH 


Surface 


7.59 


3.16 


X 


10^ 


2.43 X 10-^ 


Yes 


No 


N2H2 


Gas 


7.44 


1.00 


X 


10^ 


5.57 X 10-12 


No 


No 


CH3CN 


Surface 


7.39 


1.00 


X 


10^ 


7.14 X 10-11 


Yes 


No 


HON 


Surface 


7.37 


1.78 


X 


10^ 


1.24 X 10-^ 


Yes 


No 


CH3 


Surface 


6.96 


5.62 


X 


10^ 


4.42 X 10-12 


Yes 


No 


CH2 


Surface 


6.96 


5.62 


X 


10^ 


4.46 X 10-12 


Yes 


No 


HNC 


Surface 


6.82 


1.78 


X 


10^ 


6.08 X 10^*^ 


Yes 


No 


CH 


Surface 


6.74 


5.62 


X 


10^ 


3.99 X 10-12 


Yes 


Maybe 


C 


Surface 


6.74 


5.62 


X 


10^ 


4.03 X 10-12 


No 


No 


C2H6 


Surface 


6.63 


3.16 


X 


10^ 


1.25 X 10-^ 


No 


No 
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Table 2 — Continued 







Maximum 


Time 


Relative 


Gas-Phase 


Surface 


Species 


Phase 


Difference 


(years) 


A 1 1 

Abundance 


Detectioni 


Detection 


N 


U.1 LcXiL-Kj 


U.U ± 


■i 69 y 1 0^ 


1 QQ X 10~^2 

X . C/ 1/ A X U 


Nn 

1 >l U 


Nn 


NH 

i>l 11 


U.1 IdCC 




69 V 1 0^ 


9 09 X 1 0~^2 

Zi.UZi A XU 


1 Co 


Nn 


IN £12 


U.1 IdCC 


u.uu 


^ R9 V 1 

(J.UZi A XU 


9 HQ V 1 n~12 

Zi.yJiJ A XU 


1 


Nn 


o 




6 40 


^ 62 X 1 0^ 


1 48 X 1 0^^^ 


Nn 


Nn 

i ^ U 


C4H2 


yffi r'P 


u.ou 


^ 69 X 1 0^ 


6 67 X 1 0^1^ 

U .U 1 A XU 


Yps 

1 Co 


Nn 


C9H2 


kJ Lil IdL/C 




^ 69 X 1 0^ 


1 00 X 1 0~^2 

X . UU /\ X u 


Nn 


Nn 




kJ Lil iClL.(^ 


u .uu 


1 78 X 1 0"^ 


3 7^ X 1 0^^ 

. 1 U A X U 


1 Co 


iVldj uc 


-'--'^2^3^ 


Gas 




3 1 6 X 1 0"' 


1 01 X 10~^2 

J- .U J- /\ J- u 


Yes 

± CO 


Nn 




vjrdo 


.00 


3 1 6 X 10^ 

0. XU /N XU 


1 QO X 10~^2 

X.t7U A XU 


Nn 


Nn 




Ul LcXiL-KJ 


■i 8fi 


3 1 6 X 10^ 

(J . X U /\ X u 


3 "13 X 10~^° 

O.UO A XU 


Yp9 
1 Co 


Nn 


C3H2 


kJ U.1 Id'^C 




5 62 X 1 0^ 

. U^ /\ X u 


1 50 X 10~^ 

i .UU /\ X u 


Yes 

± CO 


Nn 

i ^ U 


OH 
vyii 


kJ Ul lCt*^C 


■i 4Q 


■i 69 X 1 0^ 


1 94 X 10~^^ 

X .^rt A X U 


Yp9 
1 Co 


Nn 


PrHn 


TTQ r*a 
kj Ul iduc 




R9 V 1 

O.UZi A XU 


R Q9 V 1 n~i2 

U . C/Zi A X u 


Nn 


Nn 


V7ll2i^ 11 


vjrclh 




1 nn V 1 

X .UU A X U 


4 1 R V 1 n~i^ 

^ . XU A XU 


Yps 

1 L.O 


Nn 

IN U 


On 

W2 


kJ Ul iclL.c 




1 78 X 1 0"^ 

X . 1 A X U 


3 97 X 1 0"^'^ 

. C/ 1 A XU 


Yps 

1 Co 


Nn 


Oo 

•^3 


vjrdo 




1 00 X 1 0^ 

X . UU /N X u 


6 71 X 1 O"^'^ 

U . 1 X A XU 


Nn 


Nn 

LWJ 




Ul IdL^c 


0. / 


IR V ID'* 

. XU A XU 


X.UO A XU 


Nn 


Nn 

IN U 


NoO 

i>i 2Vy 


Oas 

VJ CIO 


48 


1 00 X 1 0^ 

J. . UU /\ i u 


Ft 05 X 1 0~^^ 

U.UU /\ iU 


Yps 

± CO 


Nn 

i ^ U 


HnS 

ii20 


Ul lCl*^C 


41 


■i 69 X 1 0^ 

(J . U.^ A X U 


5 08 X 1 0~9 

U.UO A XU 


Yps 

1 Co 


Nn 

LWJ 


HoOo 

112^2 


Oas 


0.0 / 


3 1 6 X 10^ 


Q 99 X 1 0~^^ 


Nn 


Nn 




Vjrdb 


O.ZD 


i .UU X iU 


z . iy X iu 


Nr. 


l\Tn 
iNO 


N2O 


Surface 


3.19 


1.78 X 10^ 


1.40 X 10-12 


Yes 


No 


NO2 


Gas 


3.04 


1.00 X 10^ 


2.27 X 10-11 


No 


No 


HNO 


Gas 


2.98 


1.00 X 10^ 


2.97 X 10-1° 


Yes 


No 


HNO 


Surface 


2.76 


1.00 X 10^ 


8.76 X 10-9 


Yes 


No 


CH3CN 


Gas 


2.71 


1.00 X 10^ 


6.84 X 10-12 


Yes 


No 


C6H2 


Surface 


2.69 


3.16 X 10^ 


1.17 X 10-12 


Yes 


No 


NO2 


Surface 


2.62 


3.16 X 10^ 


1.05 X 10-12 


No 


No 


C3H 


Surface 


2.44 


1.00 X 10^ 


1.02 X 10-12 


Yes 


No 


CO2 


Surface 


2.23 


1.00 X 10^ 


2.08 X 10-1° 


Yes 


Yes 
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Table 2 — Continued 







Maximum 


Time 


Relative 


Gas-Phase 


Surface 


Species 


Phase 


Difference 


(years) 


Abundance 


Detectioni 


Detection2 


CH3CO+ 


Gas 


2.16 


1.00 X 10^ 


2.59 X 10^12 


No 


No 


HC2NC 


Gas 


2.09 


1.00 X 10^ 


1.51 X 10^12 


Yes 


No 


C3H3+ 


Gas 


2.09 


1.00 X 10^ 


2.55 X 10-12 


No 


No 


HC3N 


Gas 


2.02 


1.00 X 10^ 


2.99 X 10-11 


Yes 


No 


H2CS 


Surface 


0.141 


1.00 X 10^ 


1.93 X 10-12 


Yes 


No 


C2S 


Surface 


0.210 


5.62 X 10^ 


1.25 X 10-12 


Yes 


No 


C5H4 


Surface 


0.292 


1.00 X 10^ 


2.82 X 10-12 


No 


No 


C3H4 


Surface 


0.323 


1.00 X 10^ 


4.01 X 10-9 


No 


No 


CS 


Surface 


0.325 


1.00 X 10^ 


1.71 X 10-12 


Yes 


No 


C6H4 


Surface 


0.349 


1.78 X 10^ 


4.65 X 10-13 


No 


No 


HS2 


Gas 


0.384 


5.62 X 10^ 


6.88 X 10-12 


No 


No 


H2S2 


Gas 


0.386 


5.62 X 10^ 


6.99 X 10-12 


No 


No 


H2S 


Gas 


0.392 


3.16 X 10^ 


1.80 X 10-1° 


Yes 


No 


H3S+ 


Gas 


0.413 


5.62 X 10^ 


7.64 X 10-13 


No 


No 



Note. — We list the maximum difference factor in abundance when DR is included 
in the OSU gas-grain code relative to when it is left out. At the relevant time in the 
model evolution, we list the abundance of each species relative to total H nuclei density 
when DR is included. We also list whether there has been an interstellar detection of 
the molecule in the gas phase or on dust surfaces. 

ihttp: / / astrochymist.org/ astrochymist_ism.html 



iGibb et all (120041 ) 
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Fig. 1. — The RR rate coefficient for Mg"*" forming Mg. The dotted hne is th e rate coefficien t 
from the UMIST database and the sohd hne shows the modern calculation of lBadnelll (120061 ). 
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Fig. 2. — The DR and RR rate coefficients for C"*" forming C. The dotted hne is the RR rate 
coefficient from the UMIST database. The dashed hne shows the modern RR calculation of 



Badnelll (l2006l ) and lies almost e xactly on top o f the U MIST data. The dot-dot-dot-dashed 



line shows the DR calculation of iBadnell et al.l (120031 ). The solid line is the total RR-I-DR 
recombination rate coefficient of Badnell. 
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Fig. 3. — Same as Fig. [2] but for N"*^ forming N. 
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Fig. 4. — The abundance of surface C with respect to the total H nuclei density as a function 
of time as calculated by the OSU gas-grain code. The sohd hne is with the new RR and DR 
data included in the model and the dotted line is without the new RR and DR included. 




Fig. 5. — Same as Fig. H] but for surface CH3OH. 
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Fig. 6. — Same as Fig. H] but for surface O3. 
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Electrons 




2 4 6 8 

log time (yrs) 

Fig. 7. — Same as Fig. H] but for free electrons. 
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^ Carbon 
^ -5.0 1 ' ' ' ' ' ' ' ' ' ' ^ 




o - / 

-7.0 , , I , , , I , , , I , , , I , , ^ 

0.0 0.2 0.4 0.6 0.8 1.0 

visual extinction 

Fig. 8. — The column density of C relative to the total H nuclei column density as a function 

of visual extinction (i.e., depth into the cloud) calculated by the Meudon PDR code for a 
diffuse PDR. The solid line is with the new RR and DR included in the model and the 
dashed line is without the new RR and DR included. 
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Fig. 9.— 



Same as Fig. [H] but for a dense PDR. 



